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Abstract: We discuss the concept of discrete scale invariance and how it leads 
to complex critical exponents (or dimensions), i.e. to the log-periodic corrections to 
scaling. After their initial suggestion as formal solutions of renormalization group 
equations in the seventies, complex exponents have been studied in the eighties in 
relation to various problems of physics embedded in hierarchical systems. Only re- 
cently has it been realized that discrete scale invariance and its associated complex 
exponents may appear "spontaneously" in euclidean systems, i.e. without the need 
for a pre-existing hierarchy. Examples are diffusion-limited-aggregation clusters, rup- 
ture in heterogeneous systems, earthquakes, animals (a generalization of percolation) 
among many other systems. We review the known mechanisms for the spontaneous 
generation of discrete scale invariance and provide an extensive list of situations where 
complex exponents have been found. This is done in order to provide a basis for a 
better fundamental understanding of discrete scale invariance. The main motivation 
to study discrete scale invariance and its signatures is that it provides new insights 
in the underlying mechanisms of scale invariance. It may also be very interesting for 
prediction purposes. 



11 INTRODUCTIONI 



2 WHAT IS DISCRETE SCALE IN VARIANCE (DSlT? 





3 


WHAT ARE THE SIGNATURES OF DSI? 




3.1 Log-pcrioclicity and complex exponents 




3.2 Higher order log-periodic harmonics| 




3.3 A simpi 


e worked-out example : the Weierstrass-Mandelbrot fractal 


function 





4 


WHAT IS THE IMPORTANCE AND USEFULNESS OF DSI? 






4.1 Existence of relevant length scales 








4.2 Non-unitary field theories 










4.3 Prediction 



5 SCENARIOS LEADING TO DSI 

5.1 Built-in geometrical hierarchy 

5.1.1 Potts model on the diamond lattice 

5.1.2 Fixed points, stable phases and critical point 

5.1.3 Singularities and log periodic corrections 

5.1.4 Related examples in programming and number theory 

5.2 Diffusion in anisotropic quenched random lattices 

5.3 Cascade of ultra-violet instabilities: growth processes and rupture . . 

5.3.1 Log-periodicity in the geometrical properties 

5.3.2 Log-periodicity in time 

5.4 Cascade of structure in hydro dynaniici 

5.5 Deterministic dynamical systems 

5.5.1 Cascades of sub-harmonic bifurcations in the transition to chaos 

5.5.2 'i'wo-coupled anharmonic oscillators 

5.5.3 Near-separatrix Hamiltonian chaotic dynamics 

5.6 Animals 

5.7 Quenched disordered systems 



6 OTHER SYSTEMS 

O 'i'he bronchial tree 

6.2 Turbulence 

6.3 Titius-Bode law 

6.4 Cravitational collapse and black hole formation 

6.5 Spinodal decomposition of binary mixtures in uniform shear flow . . . 

6.6 Cosmic lacunarity 

6.7 Rate of escape from stable attractors 

O Interface crack tip stress singularity 

6.9 The Altes'wavelet and optimal time-scale product 

6.10 Orthouormal Compactly Supported Wavelets with Optimal Sobolev 
Regularity 

6.11 Eigenfunctions of the Laplace transform 

6.12 Life evolution 



rr APPLICATIONS! 



7.1 Identifying characteristic scales 

T72 'i'ime-to-tailure analysis 

7.3 Log-periodic antenna^ 

7.4 Optical waveguides 

8 OPEN PROBLEMS 

8.1 Nonlinear map and multicriticality . . . 

O IVlultilacunarity and quasi-log-periodicity 



873 Effect of disorder 



8.4 Averaging : grand canonical versus canonical 




8.5 Amplitude of log-periodicity 




8.6 Where to look for log-periodicity 




8.7 Prefered scaling ratio around 2'1 


8.8 Critical behavior and self-organized criticalit} 





1 INTRODUCTION 



During the third century BC, Euchd and his students introduced the concept of space 
dimension, which can take positive integer values equal to the number of independent 
directions. We have to wait until the second half of the nineteen century and the 
twentieth century to witness the generalization of dimensions to fractional values. 
The word "fractal" is coined by Mandelbrot to describe sets consisting of parts 
similar to the whole, and which can be described by a fractional dimension (see 
for a compilation of the most important reprints of mathemacatical works leading 
to fractals). This generalization of the notion of a dimension from integers to real 
numbers reflects the conceptual jump from translational invariance to continous scale 
invariance. 

The goal of this paper is to review the mathematical and physical meaning of 
a further generalization, wherein the dimensions or exponents are taken from the 
set of complex numbers [] (see Ref. for a recent review). We will see that this 
generalization captures the interesting and rich phenomenology of systems exhibiting 
discrete scale invariance, a weaker form of scale invariance symmetry, associated with 
log-periodic corrections to scaling. 

Before explaining what is discrete scale invariance, describing its signatures and 
importance and studying its mechanisms, let us present a brief historial perspective. 
To our knowledge, the first physical model, where sine log-periodic functions arose, 
was a model of shock waves in layered systems constructed by Zababakhin in 1965 
m . He studied the propagation of shock waves in a continuous medium with period- 
ically distributed properties and, by following the evolution of the initial impulse by 
numerical calculation, constructed a solution that also periodically reproduces itself 
on a variable scale, i.e. a log-periodic self-similar solution. A few years later, Baren- 
blatt and Zeldovich considered Zababakhin's model as a very natural development 
of self-similarity, in the following sense. They point out that the determination of 
exponents in self-similar solutions of the second type (i.e. not determined solely by 
dimension analysis) is deeply connected to the classical eigenvalue problem for linear 
differential operators. The point is that a propagating wave u{x — ci) becomes a self- 
similar homogeneous function U{X/T'^) by posing x = InX, t = InT, which writes 
again the fact that performing a scale transformation is the same as a translation in 
the logarithm. The nonlinear shock problem is thus connected to a eigenvalue prob- 
lem with Block-like wave solutions, which by a mechanism of exponentiation leads to 
the observed log-periodicity. Here, the existence of a preexisting linear periodicity is 
essential and the creation of log-periodicity can be seen as an exponentiation opera- 
tion of the discrete translational invariance into a discrete scale invariance, similarly 



to other mechanisms discussed below for discretized systems ^0 . 

Novikov has also pointed out in 1966 that structure factors in turbulence should 
contain log-periodic oscillations . Loosely speaking, if an unstable eddy in turbulent 
flow typically breaks up into two or three smaller eddies, but not into 10 or 20 eddies, 
then one can suspect the existence of a preferable scale factor, hence the log-periodic 
oscillations. The interest in log-periodic oscillations has been somewhat revived after 
the introduction of the renormalization group theory of critical phenomena. Indeed, 
the mathematical existence of such corrections has been discussed quite early in 



further generalization to the set of quaternions (the unique non-commutative generaUzation 
of complex numbers on the set of real numbers) does not bring any new structure. 



renormalization group solutions for the statistical mechanics of critical phase ||, |^. 
However, these log-periodic oscillations, which amount to consider complex critical 
exponents, were rejected for translationally invariant systems, on the (not totally 
correct |]10[) basis that a period (even in a logarithmic scale) implies the existence 



of one or several characteristic scales, which is forbidden in these ergodic systems 
in the critical regime. Complex exponents were therefore restricted to systems with 
discrete renormalization groups. In the eighties, the search for exact solution of the 
renormalization group led to the exploration of models put on hierarchical lattices, 
for which one can often obtain an exact renormalization group recursion relation. 
Then, by construction as we will show below, discrete scale invariance and complex 
exponents and their log-periodic signature appear. 

Only recently has it been realized that discrete scale invariance and its associated 
complex exponents can appear spontaneously, without the need for a pre-existing 
hierarchical structure. It is this aspect of the domain that is the most fascinating 
and on which we will spend most of our time. 



2 WHAT IS DISCRETE SCALE INVARIANCE 
(DSI)? 

Let us first recall what is the concept of (continuous) scale invariance : in a nutshell, 
it means reproducing itself on different time or space scales. More precisely, an 
observable O which depends on a "control" parameter x is scale invariant under the 
arbitrary change x ^ Ax Q if there is a number /i(A) such that 

0{x)=fiO{Xx) . (1) 

Eq.(P defines a homogeneous function and is encountered in the theory of critical 
phenomena, in turbulence, etc. Its solution is simply a power law 0{x) = Cx", with 
a = — j^l^, which can be verified directly by insertion. Power laws are the hallmark 

of scale invariance as the ratio = A" does not depend on x, i.e. the relative value 
of the observable at two different scales only depend on the ratio of the two scales 
0. This is the fundamental property that associates power laws to scale invariance, 
self-similarity and criticality 0. 

Discrete scale invariance (DSI) is a weaker kind of scale invariance according to 
which the system or the observable obeys scale invariance as defined above only for 
specific choices of A (and therefore fi), which form in general an infinite but countable 
set of values Ai, A2, ... that can be written as A„ = A". A is the fundamental scaling 
ratio. This property can be qualitatively seen to encode a lacunarity of the fractal 
structure 

^Here, we implicitely assume that a change of scale leads to a change of control parameter as in 
the renormalization group formalism. More directly, x can itself be a scale. 

•^This is only true for a function of a single parameter. Homogeneous functions of several variables 
take a more complex form than (|^). 

^Self-similarity is the same notion as scale invariance but is expressed in the geometrical domain, 
with application to fractals. 

^Criticality refers to the state of a system which has scale invariant properties. The critical state 
is usually reached by tuning a control parameter as in liquid-gas and paramagnetic-ferromagnetic 
phase transitions. Many driven extended out-of-equilibrium systems seem also to exhibit a kind of 
dynamical criticality, that has been coined "self-organized criticality" . 



Note that, since x ^ \x and 0{x) fiO{\x) is equivalent toy = \ogx y+log A 
and log 0{y) log 0{y + log A) + log ^, a scale transformation is simply a translation 
of logo; leading to a translation of O. Continuous scale invariance is thus the same as 
continuous translational invariance expressed on the logarithms of the variables. DSI 
is then seen as the restriction of the continuous translational invariance to a discrete 
translational invariance : log O is simply translated when translating y by a multiple 
of a fundamental "unit" size log A. Going from continuous scale invariance to DSI 
can thus be compared with (in logarithmic scales) going from the fluid state to the 
solid state in condensed matter physics! In other words, the symmetry group is no 
more the full set of translations but only those which are multiple of a fundamental 
discrete generator. 



3 WHAT ARE THE SIGNATURES OF DSI? 

3.1 Log-periodicity and complex exponents 

We have seen that the hallmark of scale invariance is the existence of power laws. 
The signature of DSI is the presence of power laws with comp /ex exponents a which 
manifests itself in data by log-periodic corrections to scaling. To sec this, consider 
the triadic Cantor set shown in figure 1. This fractal is built by a recursive process as 
follows. The first step consists in dividing the unit interval into three equal intervals 
of length I and in deleting the central one. In the second step, the two remaining 
intervals of length | are themselves divided into three equal intervals of length | and 
their central intervals are deleted, thus keeping 4 intervals of length |. The process 
is then iterated ad infinitum. It is usually stated that this triadic Cantor set has the 
fractal (capacity) dimension Dq — as the number of intervals grows as 2" while 
their length shrinks as 3~" at the n-th iteration. 

It is obvious to see that, by construction, this triadic Cantor set is geometrically 
identical to itself only under magnification or coarse-graining by factors Ap = 3^ which 
are arbitrary powers of 3. If you take another magnification factor, say 1.5, you will 
not be able to superimpose the magnified part on the initial Cantor set. We must 
thus conclude that the triadic Cantor set does not possess the property of continuous 
scale invariance but only that of DSI under the fundamental scaling ratio 3. 

This can be quantified as follows. Call Nx{n) the number of intervals found at the 
n-th iteration of the construction. Call x the magnification factor. The original unit 
interval corresponds to magnification 1 by definition. Obviously, when the magnifi- 
cation increases by a factor 3, the number Nx{n) increases by a factor 2 independent 
of the particular index of the iteration. The fractal dimension is defined as 

D=U„,i^«iW^^^!EW=!2I^«0.63 . (2) 

x-*oo inx x^o inx log 3 

However, the calculation of a fractal dimension usually makes use of arbitrary values 
of the magnification and not only those equal to x — only. If we increase the 
magnification continuously from say x — 3^ to x — 3^'^^, the numbers of intervals in 
all classes jump by a factor of 2 at a; = 3^, but then remains unchanged until x = 3^"'"^, 
at which point they jump again by an additional factor of 2. For 3^ < x < N^iji) 
does not change while x increases, so the measured fractal dimension D{x) ~ ^"j^'^j"^ 



decreases. The value D = 0.63 is obtained only when x is a positive or negative 
power of three. For continuous values of x one has 

NM=N^in)x''p(^-^y (3) 

where P is a function of period unity. Now, since P is a periodic function, we can 
expand it Fourier series 

P = > Cnexp\2nm- . (4) 

Plugging this expansion back into @, it appears that D is replaced by an infinity of 
complex values 

Dn = D + ni-^ . (5) 
log 3 

We now see that a proper characterization of the fractal is given by this set of com- 
plex dimensions which quantifies not only the asymptotic behaviour of the number of 
fragments at a given magnification, but also its modulations at intermediate magni- 
fications. The imaginary part of the complex dimension is directly controlled by the 
prefered ratio 3 under which the triadic Cantor set is exactly self-similar. Let us em- 
phasize that DSI refers to discreteness in terms of scales, rather than discreteness in 
space (eg like discreteness of a cubic lattice approximation to a continuous medium). 
If we keep only the first term in the Fourier series in (^) and insert in (^), we get 

iV,(n) = iViHx^fl + 2^cos(2n7ri^)^ , (6) 
V Co ln3 y 

where we have used c_i = ci to ensure that Nx{n) is real. Expression (|^) shows 
that the imaginary part of the fractal dimension translates itself into a log-periodic 
modulation decorating the leading power law behavior. Notice that the period of the 
log-periodic modulation is simply given by the logarithm of the prefered scaling ratio. 
This is a fundamental result that we will retrieve in the various examples discussed 
below. The higher harmonics are related to the higher order dimensions. 

It is in fact possible to obtain directly all these results from (|I|). Indeed, let us 
look for a solution of the form 0{x) = Cx". Reporting in (|l|), we get the equation 
1 = /iA". But 1 is nothing but e*^'^", where n is an arbitrary integer. We get then 

log/i .27m 

a = — - + r , 7 

log A log A 

which has exactly the same structure as (|^). The special case n = gives the usual 
real power law solution corresponding to fully continuous scale invariance. In contrast, 
the more general complex solution corresponds to a possible DSI with the prefered 
scaling factor A. The reason why (|T]) has solutions in terms of complex exponents 
stems from the fact that a finite rescaling has been done by the finite factor A. In 
critical phenomena presenting continuous scale invariance, (|T]) corresponds to the 
linearization, close to the fixed point, of a renormalization group equation describing 
the behavior of the observable under a rescaling by an arbitrary factor A. The power 
law solution and its exponent a must then not depend on the specific choice of A, 
especially if the rescaling is taken infinitesimal, i.e. A 1+. In the usual notation, 
if A is noted A = e""^^, this imphes that /i = e""^' and a = — ^ is independent of the 
rescaing factor £. In this case, the imaginary part in drops out. 



3.2 Higher order log-periodic harmonics 



Complex critical exponents can be derived from a discrete renormalization group. 
From it, we can obtain the relative amplitudes of the harmonics of the Fourier series 
expansion of the log-periodic function, which gives the leading correction to scaling. 
We can also discuss the effect of a first non-linear term in the flow map in the 
generation of log-periodic accumulation of singular points. 

A naive way of obtaining discrete scale covariance is to consider a discrete fractal 
for which only discrete renormalizations are allowed. Calling K the coupling (e.g. 
K = e'^l'^ for a spin model) and R the renormalization group map between two 
successive generations of the discrete fractal, one has 

f[K)=g[K)^-i{R(K)\ , (8) 

where / is the free energy per lattice site or bond, gvs,& regular part which is made 
of the free energy of the degrees of freedom summed over between two successive 
renormalizations, ^ is the ratio of the number of degrees of freedom between two 
successive renormalizations. The equation (|^ is solved recursively by 

OO 1 

f{K) = Y.-g\R^-\K)\ , (9) 

n=0 

where R^"^^ is the n*^ iterate of the renormalization transformation. For general 
ferromagnetic systems, R has two stable fixed points corresponding to = and 
K = OO fixed points, and an unstable fixed point bX, K = K^. which corresponds to 
the usual critical point. It is easy to see that / is singular at K^.. Call A the slope of 
the RG transformation at (A > 1 since the fixed point is unstable). Then the i*^ 
term in the series for the k^^ derivative of / at the fixed point will be proportional 
to i}^ I [if ^ giving rise to a divergence of the series for the k^^ derivative of / for k 
large enough, hence the singular behaviour. Now assume that j (x\K — K,}^ close 
to the critical point. Plugging this form in gives the constrain, since is a regular 
function, 

\ ra 

- = i. (10) 

This equation has an infinity of solutions given by 

mi = m + HQ , (11) 

with 

Inu ^ 27r , , 

and / is an arbitrary integer. Decorating the main algebraic behaviour, there is thus 
the possibility of terms of the form \K — Kc\"^ cos{lfl In \K — Kc\ + (pi), which are 
the looked for log-periodic corections. Such terms are not generally scale covariant. 
If continuous changes of scale were allowed, they would be incompatible with scale 
covariance, and thus forbidden. It is only because the changes of scale are discrete that 
they survive here : only when the change of scale is of the form \K — Kc\ X''\K — Kd 
do they transform multiplicatively, ensuring scale covariance. 



To find out more about these correction terms, it is very convenient to use tfie 
Mellin transform as done in 0. Introduce therefore for any function / the Melhn 



transform _ 

coo 

5-1 



f{s) = / x^-'f{x)dx . (13) 





We then consider the free energy in (P) in the hnear approximation, replacing i?^") (K) 
by Kc+X'^{K—Kc). Setting x = K—K^ we have then, by applying this transformation 
to both sides of (^ 

We then reconstruct the original function by taking the inverse transform 

1 rc+ioo „ 

fix) = — / f{s)x-'ds . (15) 

ZlTT Jc-ioo 

The usefulness of the Mellin transform is that the power law behavior springs out 
immediately from the poles of /(s), using Cauchy's theorem. For a general statistical 
mechanics model, g being the regular part of the free energy has generally the form 
of the logarithm of a polynomial in x. Factorizing the polynomial, we do not lose 
generality by considering g given by 

g{x) = ln(l + x) , (16) 

for which 

9is) = — ^ . (17) 
s sm SIT 

In inverting the Mellin transform, we have two types of poles. The poles of g occur 
for s = —n, n > 0, and contribute only to the regular part of /, as expected since 
g is a regular contribution. The poles of the prefactor in (|T4|), which stem from the 
infinite sum over successive embeddings of scales, occur at 

s = -mi , (18) 

with mi as in ([ll|. They correspond exactly to the singular contributions discussed 
above. Their amplitude is obtained by applying Cauchy's theorem and is of the order 
of 

1 1 



m + HQ siniT {m + HQ) ' ^ ^ 

which behaves at large I as e"'^'^. Hence the amplitude of the log-periodic corrections 
decays exponentially fast as a function of the order I of the harmonics. This explains 
the fact that mostly the first harmonic accounts for experimental data. 

Of course the previous computation suffers from the linear approximation, which 
becomes deeply incorrect as n gets large in (|^), hence in the region determining the 
singularity. As discussed in |Q, the crucial property missed by the linear approxima- 



tion is that / is analytic only in a sector | ar (7x| < 9 while we treated it as analytic in 
the cut plane \argx\ < n. The true asymptotic decay of the amplitudes of successive 
log-periodic harmonics is therefore slower than initially found, and goes as e~'^^. The 
angle 6 depends specifically on the flow map of the discrete renormalization group 
p3| and is generally of order 1. Assuming this formula e~'^^ of the decay still gives 



the right order of magnitude for / = 1, we see that the amphtude of the first mode 
can be much larger than found above in the hnear approximation. 

The hnear approximation misses more interesting features, discussed in [ |136[j . 
Beyond the hnear approximation, the renormahzation transformation R can be ex- 
panded as a polynomial in x. Let us keep the first two terms only, so R{x) = \x~vx^ 
with z/ > 0. Such a non monotonic RG transformation occurs for instance in frus- 
trated models with both ferromagnetic and antiferromagnetic interactions. An in- 
teresting consequence is that this transformation has another unstable fixed point at 
Xc = The function / is singular at all the pre-images of this other fixed point, 
if the absolute value of the Lyapunov exponent of R at this fixed point is greater 

2 

than one namely if 1 < A < 1 + | or A < A_ or A > A+ where A± = — ^ ^ — — — . 
This is due to the fact that these preimages all go to the fixed point after a finite 
number of renormalizations. Moreover these preimages accumulate at the original 
unstable fixed point x = (figure 1), and their accumulation becomes geometrical 
very close to x = 0. When crossing such a preimage, there is a singularity in /, 
usually manifested as a kink in the curve. This sequence of kinks close to x = can 
be fitted well with the first harmonic of log-periodic oscillations discussed previously. 
But the discussion suggests that there is more to be seen in these oscillations - the 
geometrical accumulation of critical points towards the "main" critical point. 



3.3 A simple worked-out example : the Weierstrass-Mandelbrot 
fractal function 

Mandelbrot proposed the following extension of the Weierstrass function (in order to 
cure the disavantage that the Weiestrass function has a larger scale) 

+00 g«7"t 

^(^) = J2 {2-D)n ^"^^ ' ^^^^ 1 < D < 2, 7 > 1 and 0„ arbitrary . (20) 

D is the Hausdorf-Besicovich fractal dimension of the graph of W{t) (i.e. the graphs 
of the real part or imaginary parts of W . We follow |12| to present some important 
properties that ihustrate the previous analysis. 

Consider the case where the phases 0^ = /in. Then, W{t) obeys the renormahza- 
tion group equation 

W{-it) = e-'^ 7^-^ Wit) . (21) 

This implies that the whole function W can be reconstructed from its value in the 
range to <t < ■jtQ. 

Note that the equation (|21|) does not imply that W{t) is a fractal, because it is 
satisfied by smooth functions of the form 

fruit) = t^'"" expH(/i + 2nm)^^] , (22) 

ln7 

where m is an integer. Wit) can thus be constructed from an infinite series of sum 
terms. The Poisson formula 



+00 +00 

E /(^) = E / fit) e'^^'* dt , (23) 



also us to obtain this new series as 

W{t)= E / ■ (24) 

After integrating by part, the terms in this series can be expressed as F functions: 
iy(t) = e^f ^ e-i^i°^ £ Ut) e— ^/"-^r(ZJ-2 + /|^'''" ) . (25) 

Apphcation of Stirhng's formula shows that the terms decay as exp(— 2m7r^/ ln7) 
as m 00 and as |m|~*^2~^) exp(2cm) (where c is a constant) as m — —00, thus 
ensuring the convergence of the series. The convergence is fastest when 7 — 1 in 
constrast to ( pOD whose convergence is fastest when 7 00. 
Consider the case /i = vr leading to 

+00 r -i\n 

A{t) = Imaginary part W{t)\^=^ = ^ -^^^ sin(7"t) . (26) 

In this case, the term m = — 1 in (|25|) dominates and gives the following large scale 
behavior 

A(t) ^ -t^-^ Imaginary part fe^'^^+T^ r(D - 2 - It^)] . (27) 

In 7 V ln7 / 

If 7r/ln7 is sufficiently large, this can be simplified by using Stirling's formula, to 
give the trend 

A{t) ^ (i^)-^ t-^ sin (7^ In f + ^) , with to = ^ . (28) 

V m7 vr Vm7 to 4/ em 7 

This is a log-periodic function with a prefered scaling ratio 7. The higher order terms 
in (^) produces the non-differentiable fractal structure of the Weierstrass-Mandelbrot 
function. 

This calculation makes explicit the relationship between the main log-periodic 
dependence and the higher order fractal structure. 

4 WHAT IS THE IMPORTANCE AND USEFUL- 
NESS OF DSI? 

4.1 Existence of relevant length scales 

Suppose that a given analysis of some data shows log-periodic structures. What can 
we get out it them? First, as we have seen, the period in log-scale of the log-periodicity 
is directly related to the existence of a prefered scaling ratio. Thus, log-periodicity 
must immediatly be seen and interpreted as the existence of a set of prefered charac- 
teristic scales forming all together a geometrical series A"^, A"^"*"^, A, A^, A", 
The existence of such prefered scales appears in contradiction with the notion that 
a critical system, exhibiting scale invariance has an infinite correlation length, hence 



only the microscopic ultraviolet cut-off and the large scale infra-red cut-off (for in- 
stance the size of the system) appear as distinguishable length scales. This recovers 
the fact that DSI is a property different from continuous scale invariance. In fact, 



it can be shown [|T4| that exponents are real if the renormalization group is a gra- 
dient flow, a rather common situation for systems at thermal equilibrium, but as 
we will see, not the only one by far. Examples when this is not the case can be 
found especially in random systems, out-of-equilibrium situations and irreversible 
growth problems. In addition to the existence of a single prefered scaling ratio and 
its associated log-periodicity discussed above, there can be several prefered ratios 
corresponding to several log-periodicities that are superimposed. This can lead to a 
richer behavior such as log-quasi-periodicity. Quasiperiodicity has been suggested to 
describe the scaling properties of diffusion-limited-aggregation clusters 

Log-periodic structures in the data indicate that the system and/or the underlying 
physical mechanisms have characteristic length scales. This is extremely interesting as 
this provides important constraints on the underlying physics. Indeed, simple power 
law behaviors are found everywhere, as seen from the explosion of the concepts of 



fractals, criticality and self-organized-criticality |Tl[]. For instance, the power law 
distribution of earthquake energies which is known as the Gutenberg-Richter law 
can be obtained by many different mechanisms and a variety of models and is thus 
extremely limited in constraining the underlying physics. Its usefulness as a modelling 
constraint is even doubtful, in contradiction with the common belief held by physicists 
on the importance of this power law. In contrast, the presence of log-periodic features 
would teach us that important physical structures, that would be hidden in the fully 
scale invariant description, existed. 

4.2 Non-unitary field theories 

In a more theoretical vein, we must notice that complex exponents do not appear 
in the canonical exactly solved models of critical phenomena like the square lattice 
Ising model or Bose Einstein condensation. This is because such models satisfy some 
sort of unitarity. From conformal invariance |T^, it is known that the exponents of 
two dimensional critical models can be measured as amplitudes of the correlation 
lengths in a strip geometry. Since the Ising model transfer matrix can be written in 
a form which is symmetric, all its eigenvalues are real, therefore all its exponents are 
real. The other standard example where exponents can be computed is e expansion. 
However in that context, there is an attitude, inherited from particle physics, to 
think mostly of Minkowski field theories. For instance in axiomatic field theory, 
Euclidian field theories are defined mostly as analytic continuations of Minkowski 



field theories. Now complex exponents, as we have argued |]10[, make perfect sense 
for Euclidian field theories, but lead to totally ill-behaved Minkowski field theories, 
with exponentially diverging correlation functions. An approach based on any sort 
of equivalence between the two points of view is bound to discard complex exponents 
(as well say as complex masses). The complex exponents can thus be viewed as 
resulting from the breaking of equivalence (or symmetry under Wick rotation) of the 
Euclidian and Minkowski field theories. As it is now undertood that quantum field 



theories are only effective theories that are essentially critical [|T^, could there be 

^The microscopic cut-off is the Planck scale ~ 10~^^m while the macroscopic cut-offs (or cor- 
relation lengths) corresponding to the observed particle masses such as for the electron are of the 



a relation between the spectrum of observed particle masses and the characteristic 
scales appearing in DSI and its variants and generalizations? 



4.3 Prediction 



Lastly, it is important to stress the practical consequence of log-periodic structures. 
For prediction purposes, it is much more constrained and thus reliable to fit a part of 
an oscillating data than a simple power law which can be quite degenerate especially 
in the presence of noise. This remark has been used and is vigorously investigated in 
several applied domains, such as earthquakes []TH|, |2T| , rupture prediction 

and financial crashes p3|, p5| . 



5 SCENARIOS LEADING TO DSI 

After the rather abstract description of DSI given above, we now discuss the physical 
mechanisms that may be found at its origin. It turns out that there is not a unique 
cause but several mechanisms may lead to DSI. Since DSI is a partial breaking of a 
continuous symmetry, this is hardly surprising as there are many ways to break down 
a symmetry. We describe the mechanisms that have been studied and are still under 
investigation. The list of mechanisms is by no mean exhaustive and other mechanisms 
may exist. We have however tried to present a rather complete introduction to the 
subject. 

It is essential to notice that all the mechanisms described below involve the exis- 
tence of a characteristic scale (an upper and/or lower cut-off) from which the DSI can 
develop and cascade. In fact, for characteristic length scales forming a geometrical 
series to be present, it is unavoidable that they "nucleate" from either a large size 
or a small mesh. This remark has the following important consequences : even if the 
mathematical solution of a given problem contains in principle complex exponents, 
if there are no such cut-off scales to which the solution can "couple" to, then the 
log-periodicity will be absent in the physical realization of the problem. An example 
of this phenomenon is provided by the interface-crack stress singularity described 
below. 



5.1 Built-in geometrical hierarchy 



The most obvious situation occurs when some physical system is put on a pre-existing 
discrete hierarchical system, such as the Bethe lattice, or a fractal tree. Since the 
hierarchical system contains by construction a discrete hierarchy of scales occurring 
according to a geometrical series, one expects and does find complex exponents and 
their associated log-periodic structures. Examples are fractal dimensions of Cantor 
sets p6|, |2^, percolation ||29|, ultrametric structures wave propagation in 



fractal systems |^, magnetic and resistive effects on a system of wires connected 
along the Sierpinski gasket [p^ , Ising and Potts models ||3^, fiber bundle rupture 
p5| , |20[| , sandpiles [^. Quasi-periodic and aperiodic structures can also often be 



captured by a discrete renormalization group and can be expected to lead to discrete 



order of 10 ^^m. This is a situation where the correlation length is thus lO'^^ times larger than the 
"lattice" size, very close indeed to criticality! 



scale invariance and log-periodicity. For instance, the quantum XY spin chain with 
quasi-periodic two-valued exchange couplings [^] has its zero-field specific heat and 



susceptibility behaving as power laws of the temperature, modulated by log-periodic 
modulations. The zero-temperature magnetization is a Cantor function of the applied 
field and the density of state is also a power law of the distance to gaps with log- 
periodic modulations. Similarly, log-periodic corrections to scaling of the amplitude of 
the surface magnetization have been found for aperiodic modulations of the coupling 
in Ising quantum chains [Q. Analytical and numerical calculations associated with 
succesive hierarchical approximations to multiscale fractal energy spectra show that, 
in a certain range of temperatures, the specific heat displays log-periodic oscillations 
as a function of the temperature ||3^. The scattering cross-section of a radiation 
scattered off a DSI fractal structure also exhibits log-periodic oscillations as a function 



of the scattering wavenumber parameter [40 



Multifractal models of geophysical fields have used discrete scale ratios to con- 
struct isotropic and anisotropic cloud and rain structures El], H2 . 



5.1.1 Potts model on the diamond lattice 

Let us now give some details to see more clearly how physics on hierarchical systems 
leads to log-periodicity. As a canonical example, we treat the Potts model on 
the diamond lattice [Q. This lattice is obtained by starting with a bond at magni- 
fication 1, replacing this bond by four bonds arranged in the shape of a diamond at 
magnification 2, and so on, as illustrated figure 2. At a given magnification 2^, one 
sees 4P bonds, and thus |(2 + 4^) sites. 

The spins o"j are located at the vertices of the diamond fractal. In the same way 
that the lattice appears different at different scales from a geometrical point of view, 
one sees a different number of spins at different scales, and they will turn out to 
interact in a scale dependent way. For a given magnification x = 2^, the spins we can 
see are coupled with an interaction energy 

E = -JJ2 ^i^^^j)^ (29) 

<ij> 

where J is the coupling strength, the sum is taken over nearest neighbors and the delta 
function equals one if arguments are equal, zero otherwise The system is assumed at 
thermal equilibrium, and the spin configurations evolve randomly in time and space 
in response to thermal fiuctuations with a probability proportional to the Boltzmann 
factor e~^^ , where (3 is the inverse of the temperature. The partition function Z at 
a given magnification x = 2^ is 

where the sum is taken over all possible spin configurations which can be seen at that 
scale. We do not compute Zp completely, but first perform a partial summation over 
the spins seen at one scale and which are coupled only to two other spins. This is 
how, in this particular example, one can carry out the program of the renormalization 
group by solving a succession of problems at different scales. Let us isolate a particular 
diamond, call ai, ct2 the spins at the extremities and Si, S2 the spins in between as in 
figure 2. The contribution of this diamond to e~^^ is 

J^S{<Ti ,si)+5{a2 ,si)+5{<Ti ,S2)+S{a2 ,S2) 



where we have defined K 



Since si,S2 enter only in this particular product, 



we can perform summation over them first when we compute Zp. The final result 
depends on whether ai and cr2 are equal or different: 



CTl,Sl) + 5{cT2,Sl)+S(ai,S2) + 5(cT2,S2) 



Sl,S2 



{2K + Q 



2)^ 



0-1 7^ ^2 



(30) 
(31) 



so we can write 



o-i,Si)+(5((T2,Sl)+5(cri,S2)+(5((T2,S2) 



Sl,S2 



^ ^ ^ [ \{2K + Q~2y 
{2K + Q -2fK'^^''^^"^^ , 



(32) 
5(0-1,0-2) 



where we used the identity 
and we set 

K' = 



2K + Q-2 



(33) 



(34) 



If we perform this partial resummation in each of the four diamonds, we obtain exactly 
the system at a lower magnification x = 2^~^. We see therefore that the interaction of 
spins tranforms very simply when the lattice is magnified : at any scale, only nearest 
neighbor spins are coupled, with a scale dependent coupling determined recursively 
through the renormalization group map 



K 



v-i 



'Kl + Q-V 
2Kp + Q-2^ 



(35) 



The spins which are "integrated out" by going from one magnification to the next 
simply contribute an overall numerical factor to the partition function, which is equal 
to the factor {2K + Q — 2)^ per edge of (0). Indeed, integrating out the spins si 
and S2 leaves only ai and a2 whose interaction weig ht is by definition K'^(''^'''^\ if 
K' represents the effective interaction weight at this lower magnification 2^^^. The 
additional numerical factor shows that the partition function is not exactly invariant 
with the rescaling but transforms according to 



Zp{K) = Zp^i[4>{K)]{2K + Q-2 



>2.4P 



(36) 



since there are 4^ bonds at magnification 2^. Now the free energy, which is defined 
as the logarithm of the partition function per bond, reads 



fpiK) 

From (|36D, we deduce the following 



4P+1 



In Zp{K) 



aK) = g{K) + \u_,{K') , 



(37) 



where 

giK) = ^lni2K + Q-2) . (38) 

For an infinite fractal, tlie free energy for some microscopic coupling K satisfies 
therefore 

fiK) = g{K) + -f{K') , (39) 

where = 4. This explicit calculation makes clear the origin of the scaling for the free 
energy : the interaction weights remain of the same functional form at each (discrete) 
level of magnification, up to a multiplicative factor which accounts for the degrees of 
freedom "left-over" when integrating from one magnification to the next. This is the 
physical origin of the function g in (|39|). 



5.1.2 Fixed points, stable phases and critical point 

Consider the map K' = (j){K) (|35|) . It exhibits three fixed points (defined by 
K' = K = (piK)) located at = 1,K = oo,K = where is easily deter- 
mined numerically, for instance ^ 3.38 for Q = 2, K^. ^ 2.62 for Q = 1. That 
K = 1 and K = oo are fixed points is obvious. The former corresponds to totally 
uncoupled spins, the latter to spins which are forced to have the same value. In 
both cases, the dynamics disappears completely, and one gets back to a purely ge- 
ometrical problem. Observe that these two fixed points are attractive. This means 
that if we start with some coupling say K > Kc deep down in the system, that is 
for very large magnifications, when one diminishes the magnification to look at the 
system at macroscopic scales, spins appear almost always parallel, and therefore are 
more and more correlated as one reduces magnification. Similarly if we start with 
K < Kc spins are less and less correlated as one reduces magnification. The condition 
K > Kc together with the definition K = e^"^ implies /3 > /3c, i.e. corresponds to 
the low-temperature regime dominated by the energy. The physical meaning of the 
attraction of the renormalization group flow to the fixed point K = oo, i.e. zero 
temperature, means that the macroscopic state of the spins is ferromagnetic with a 
macroscopic organization where a majority of spins have the same value. Similarly, 
the condition K < Kc implies P < Pc, i-e. corresponds to the high-temperature 
regime dominated by the entropy or thermal agitation. The physical meaning of 
the attraction of the renormalization group flow to the fixed point K = 0, i.e. infi- 
nite temperature, means that the macroscopic state is completely random with zero 
macroscopic magnetization. 

The intermediate fixed point Kc, which in contrast is repulsive, plays a completely 
different and very special role. It does not describe a stable thermodynamic phase 
but rather the transition from one phase to another. The repulsive nature of the 
renormalization group map flow means that this transition occurs for a very special 
value of the control parameter (the temperature or the coupling weight K = Kc). 
Indeed, if we have spins interacting with a coupling strength right at Kc at microscopic 
scales, then even by reducing the magnification we still see spins interacting with a 
coupling strength right at Kc\ This is also a point where spins must have an infinite 
correlation length (otherwise it would decrease to zero as magnification is reduced, 
corresponding to a different effective interaction): by definition it is a critical point. 



Close to Kc we can linearize the renormalization group transformation 



K'-Kc^ X{K - Kc) , (40) 

where X = -^{k^ > 1. For couplings close enough to the critical point, we now see 
that as we increase magnification, the change in coupling becomes also very simple ; 
only, it is not the coupling that gets renormalized by a multiplicative factor, but the 
distance to Kc. 

The equation (|39|) together with ( ^0|) provides an explicit realization of the pos- 
tulated functional form (up to the non-singular term g), where the coupling 
parameter K (in fact K — Kc) plays the role of the control parameter x. 



5.1.3 Singularities and log periodic corrections 

The renormalization group equations (^) and (^) can be solved for the free energy 
by 

OO -| 

fiK) = T.-^9[<P^''\K)] , (41) 

n=0 

where (^'•"^ is the n*^ iterate of the transformation (J) (eg 0*^^^(x) = 0[(/)(x)]). Now it 
is easy to show [ p.36| | that the sum (^) is singular at i^' = Kc- This stems from 
the fact that Kc is an unstable fixed point, so the derivative of at i^c is A > 1. 
Therefore if we consider the k^^ derivative of / in (|4l|) it is determined by a series 

which is greater than 1 for k large enough, so 
this series diverges. In other words high enough derivatives of / are infinite at Kc- 
Very generally, this implies that close to Kc one has 

f{K)^{K-Kc)'^ , (42) 

where m is called a critical exponent. For instance if < m < 1, the derivative 
of / diverges at the critical point. Plugging this back in (|39D, we see that, since g 
is regular at Kc as can be checked easily from (|38D, we can substitute it in (^) and 
recover the leading critical behavior and derive m solely from the equation {K — 
Kc)'^ = j^[X{K — Kc)]"^ involving only the singular part, with the fiow map which has 
been linearized in the vicinity of the critical point. Therefore, the exponent satisfies 
A™ = /i, an equation that we have already encountered and whose general solution is 
given by 

m A m A 

To get expression (|43|), we have again used the identity e*^'"" = 1. We see that 
because there is discrete scale invariance (namely (^) holds which relates the free 
energy only at two different scales in the ratio 2), nothing forces m to actually be 
a real number. In complete analogy with the case of complex fractal dimension, a 
critical phenomenon on a fractal exhibits complex critical exponents. Of course 
/ is real, so the most general form of / close to the critical point should be 



f{K) ^{K- Kc)"" jao + E a„ cospvml^ \n{K - Kc) + | 

I n>0 J 



(44) 



where 

In u ^ 1 , , 



hence exhibiting the log-periodic corrections. Derrida et aL have studied this 
example more quantitatively and find that the amplitude of the log-periodic oscil- 
lations are of the order of 10"*^ times less that the leading behavior. This is thus a 
small effect. In contrast, the examples below exhibit a much stronger amplitude of 
the log-periodic corrections to scaling, that can reach 10% or more. 



5.1.4 Related examples in programming and number theory 

Log-periodicity, many of which are of a fractal nature, are found in the solutions of 



algorithms based on a recursive divide- and- conquer strategy |^4[ such as heapsort, 
mergesort, Karatsuba's multiprecision multiplication, discrete Fourier transform, bi- 
nomial queues, sorting networks, etc. For instance, it is well-known that the worst 
time cost measured in the number of comparisons that are required for sorting n 
elements by the MergeSort procedure is given by n log2 n to leading order. It is less 



known that the first subleading term is nP(log2n), where P is periodic K 



Reducing a problem to number theory is like striping it down to its sheer fun- 
damentals. In this vein, arithmetic functions related to the number representation 
systems exhibit various log-periodicities. For instance, the total number of ones in 
the binary representations of the first n integers is log2 n + nF {\og2 n), where F is 



a fractal function, continuous, periodic and nowhere differentiable P5| . 

The statistical distribution of energy level spacings in two-dimensional harmonic 
oscillators with irrational frequency ratio uji/uj2 also exhibit a discrete hierarchical 
structure deeply connected to the properties of irrational numbers. This problem 
arises in the study of the correspondence between classical and quantum dynamics in 
Hamitonian systems [|^. For integrable classical systems, the generic level spacing 



distribution is the Poisson's law of exponential decay with a maximum at zero, corre- 
sponding to level clustering. The two-dimensional harmonic oscillator problem, pos- 
sibly the simplest integrable system, does not follow this generic rule. The spectrum 
is found to exhit a discrete and rigid structure controlled by the continued fraction 



expansion of the irrational number uJi/iU2 0|- Following the Einstein-Brillouin-Keller 



(EBK) semi-classical quantization rule, which for harmonic oscillators give the exact 
quantum mechanical results, the energy levels are 

hi 1 

Emi,m2 = TT- li^l + o)^l + + -)^2] ■ (46) 

In units of h/2n, we see that the difference between two energy levels is {■mi—m[)uji + 
(m2 — 7712)002- For UJ1/UJ2 irrational, this is never zero. Consider first the simplest and 
best irrational, the golden mean ai = uji/uj2 = (\/5 — l)/2 which has the continued 
fraction expansion ai = [1, 1, 1, ...]. The closest difference are obtained for mi — m[ 
and mg— m2 being the successive Fibonacci numbers (Fq = Fi = 1, Fn+i = F/v+F„_i) 
which are such that Fn/Fn+i are the best rational approximations to ai. Due to their 
property and that of the golden mean, we have — ctiF„ = (— cri)"+^. Thus all 
level spacings are integer powers of cxi. The entire distribution of level spacings AE 
satisfies 

P,,(AF) ~ 5(AE -O , for n = 1,2,... (47) 

We recognize a special case of a discrete scale invariant distribution with prefered 
scaling ratio Ui. 



This results holds for other irrational numbers. Consider uji/uj2 = o"2 = "\/2 — 1 = 
[2, 2, 2, ....]. Now, P„^{^E) not only peaks at /S.E = \GnCr2 — Gn-i correspond- 
ing to the continued fraction approximations Gn-i/Gn (with Gq = Gi = 1 and 
Gn+i = 2Gn + Gn-i), it also shows peaks that correspond to the so-called interme- 



diate fractions {Gn + G'„,-i)/(G„+i + Gn) [ill. These intermediate fractions are the 



second best approximations. The peak corresponding to the intermediate fraction 
{Gn + Gn-i)/{Gn+i + Gn) has the same height as the one corresponding to Gn-i/Gn, 
while these main peaks continue to follow the inverse- square law. 

For a general ratio lji/u2 = [01,02,03, ....,a„, ...], it is best approximated by the 
successive fractions p/q = [ai, 02, as, aj obtained by successive truncations of the 
continued fraction expansion. In addition, when a„ > 1, there are a„ — 1 intermedi- 
ate fractions (j9„_i + kpn)/{qn-i + kqn) with k = 1,2, ...,an — 1, between Pn-i/q-n-i 
and Pnlq-n which also provide a good approximation to uji/uj2. There is a one-to-one 
correspondence between the allowed nearest neighbor level spacings and the contin- 



ued and intermediate fraction approximations of the frequency ratio uji/uj2 [^. In 
the case where uji/uj2 belongs to the class of relatively simple irrationals known as 
quadratic numbers, for which the continued fraction expansion is periodic (which in- 
cludes o"i and (J2), the entire distribution P^^^i^^^i^^E) is self-similar and log-periodic: 
/S.E'^ P^-^^/^^{AE) is periodic in a log-log scale with period given by the logarithm of 
the prefered scaling ratio c<Ji/u;2. This log-periodicity is obviously related to the fact 
that the continued fraction expansion of a quadratic irrational number is periodic. 
Since going from one level to the next in the continued fraction expansion is a mul- 
tiplicative operation, the periodicity in the continued fraction expansion is reflected 
in a log-periodicity in the number properties. 

The exact log-periodicity is lost for more general irrational numbers while the 
general scaling inverse square law seems to persist. In addition, local log-periodicity 
can be observed resulting from local periodicity of the continued fraction expansion 
of the irrational number. In sum, log-periodicity arises in this problem because one 
looks for an extremal property, namely the successive spacings, in the presence of a 
periodic continued fraction expansion. 

5.2 Diffusion in anisotropic quenched random lattices 

In this scenario, the DSI hierarchy is constructed dynamically in a random walk 
process due to intermittent encounters with slow regions [Q. Consider a random 
walker jumping from site to site. Bonds between sites are of two types : (i) directed 
ones on which the walker surely goes from his site to the next on his right ( "diode" 
situation) ; (ii) two-way bonds characterized by a rate u (resp. v) to jump to the 
neighboring site on his right (resp. left). The fraction of two-way bonds is 1 — p 
and the fraction of directed bonds is p. We construct a frozen random lattice by 
choosing a given configuration of randomly distributed mixtures of the two bond 
species according to their respective average concentration p and 1 — p. The exact 
solution of this problem has been given in |ED[ and shows very clearly nice log-periodic 



oscillations in the dependence of {x ) as a function of time, as seen in figure 3. 



We now present a simple scaling argument which recovers the exact results. 
To do so, we assume - << 1 — »<< 1. We are thus in a situation where most bonds 
are directed and dilute clusters of two-way bonds are present. In addition, the two 
bonds are strongly impeding the progress of the walker as the forward rate to the 



right is much smaller than the backward rate to the left. 

In this situation, the random walker progresses at constant velocity to the right as 
long as it encounters only diode bonds and gets partially trapped when it encounters 
two-way bonds. To see how DSI is spontaneously generated, we estimate the typical 
number of jumps needed for the random walker to pass k adjacent two-way bonds, 

i.e. a connected cluster of k two-way bonds. In the limit f << 1, n ~ (^) [|. 
Using the fact that the average separation between fc-tuples of two-way bonds is 
approximately (1 — if ^ << (1 — p), the typical number of jumps needed for 
the random walker to go beyond the first fc-tuple of consecutive two-way values is 
completely dominated by r^. One thus expects the rate as a function of the number of 
jumps to exhibit local minima at r = r^: these are the jump-numbers time scales. The 
second part of the scaling argument consists in recognizing that the random walker has 
to cover a typical distance from the origin to encounter the first k-tuple of consecutive 
'two-way' values, of the order of (1 —p)~^. In the limit where ^ << (1 — p) « 1, 
we can thus write the approximate renormalization group equation 

x{T)^{l~p)x{\r)+g{r) , (48) 

where we have set A = ^ and (^(r) is some regular function taking into account various 
local effects that correct the main scaling. Notice the similarity with (|39D . Because 
this renormalization group equation can be written only at scales which are powers 
of A, we are back to the situation discussed before. We see in particular that x ^ 
with 

„=!2l(i^ + i^. (49) 

\.og[u/v) iogiu/v) 

This result for the exponent v turns out to be exact. Furthermore, the range of 
parameters over which this holds is much larger than suggested by this intuitive 
argument. More precisely, as soon as - < (1— p), one finds x{t) ~ P(i^gf^), where 



log(l-p) 



iog(w/D) ^^"^ -P is a periodic function of unit period. This prediction is remarkably 



well-confirmed by numerical simulations and recovers the exact calculation of . 
This is shown in figure 3. 

These log-periodic oscillations are not artifacts of ID-systems as random walks 
with a fixed bias direction on randomly diluted 3D-cubic lattices far above the per- 
colation threshold have been found to exhibit log-periodic oscillations in the effective 



exponent versus time ||5^. These log-periodic oscillations stem from an Arrhenius 
exponentiation of a discrete spectrum of trapping well depths, leading to a discrete 
hierarchy of time scales. This mechanism is very similar to the previous one. Similar 
results have been found also for two-dimensional systems under the same conditions 
57| . Julia Draeger in Hamburg has also found log-periodic oscillations above the two- 



dimen- sional percolation tthreshold with topological bias in diffusion, emphasized by 
Havlin and Bunde a decade ago. The topological bias has not a fixed direction like 
an electric field, but points along the local direction of the shortest distance (more 
appropriate for hydrodynamic flow). 

This mechanism for generating log-periodic oscillations makes use of an interplay 
between dynamics and quenched randomness leading to a regime where the dynamics 

^and not Tk ~ fc^. This is due to the fact that the walker goes back many times before escaping 
from the cluster of size k. 



is highly intermittent. The presence of the discrete lattice and the mesh size is es- 
sential. This can be generalized : when the property of a system depends on extreme 
fluctuations which have a size or amplitude equal to a multiple of a characteristic 
mesh size, log-periodic oscillations can be expected. This mechanism thus applies 
also to the Lifshitz singularities in random systems : exact calculations in random 
harmonic one- dimensional chains exhibit a (log-)periodic factor in the Lifshitz essen- 
tial singularity occuring near the van Hove band-edge in the density of state fSTJ. 



Intuitively, in the language of harmonic chains [^, this essential singularity stems 
from the fact that the existence of a mode with a frequency close to the maximal 
frequency demands a large region containing light masses only. Similar intermittent 
amplification processes have also been studied in |53|. Log-periodicity found in 



the solutions of boolean delay equations p3, K3| stems from a similar mechanism. 



Let us also mention Ref. ||58| who report log-periodic oscillations in the survival 
probability for particles at the interface between an accelerator mode and the chaotic 
sea in the kicked rotor probelm. The log-periodicity seems to be reinforced by quan- 
tum effects. Its origin could be the existence of a regular time scale, the kicking time, 
which plays the role of a ultraviolet cut-off and the anomalous trapping/ acceleration 
leads to an exponentiation of this arithmetic (i.e. additive) discrete time scales. Thus 
discrete scale invariance comes here from the same type of exponentiation of a discrete 
translational invariance, as before. 

5.3 Cascade of ultra-violet instabilities : growth processes 
and rupture 

5.3.1 Log-periodicity in the geometrical properties 

Numerical analysis of large diffusion-limited-aggregates have uncovered a discrete 
scaling invariance in their inner structure, which can be quantified by the introduc- 
tion of a set of complex fractal dimensions . The values of the complex fractal 
dimensions can be predicted quantitatively from a renormalization group approach 
using the quasi-periodic mapping found in . 

A theoretical investigation of a simplified model of DLA, the needle problem, 
which is also of direct application to crack growths has been done to identify the 
underlying physical mechanism |6^. Based on perturbative analysis and some exact 



results from the hodograph method in the 2D conformal plane, we find that the two 
basic ingredients leading to DSI are the short wavelength MuUins-Sekerka instability 
f\ and the strong screening of competing needles. The basic simple picture that 
emerges is that non-linear interactions between the unstable modes of the set of 
needles lead to a succession of period doubling, the next sub-harmonic catching up 
and eventually screening the leading unstable mode. The succession of these period 
doubling explains the existence of discrete scale invariance in these systems. We thus 
think that short wavelength instabilities of the Mullins-Sekerka type supplemented 
by a strong screening effect provides a general scenario for the spontaneous formation 
of log-periodic structures. This scenario provides, in addition, an explanation for the 
observation of a prefered scaling ratio close to 2. 

^The Mullins-Sekerka instability is nothing other than the "lightning rod effect" well-known in 
electrostatics, according to which large curvature concentrate the gradient of the potential field. 
Here, the growth velocity is proportional to the gradient of the concentration field. 



Numerical simulations of the needle problem, using various growth rules (DLA, 
angle screening, r^-model, crack approximation) on systems containing up to 5000 
needles confirm clearly the proposed scenario, as shown in figure 4. The density of 
needles as a function of the distance to the base presents clear evidence of log-periodic 
modulations of the leading algebraic decay. Geological data on joints competing in 
their growth in a similar fashion also exhibit approximately the log-periodic structure 
|]60| , |6l| . |]62| , |63| present further data on joints which, in our eyes, exhibit clearly 
log-periodicity, even if the authors were not aware of the concept. Various previous 
investigation of the growth of arrays of cracks have shown the log-periodic structures, 
even if the authors neither point it out nor explained the mechanism . 

What we learn by comparing these different systems, with various growth rules, is 
that the spontaneous formation of DSI seems robust with respect to significant modi- 
fications. The improvement of our understanding of DLA resides on the identification 
of a spontaneous generation of an approximately discrete cascade of Mullins-Sekerka 
instabilities from small scales to large scales. This discreteness results from a cascade 
of mode selections by a nonlinear nonperturbative coupling between modes of growth 



60 1 . Let us mention that, in the early eighties, Sadovskiy et al. have argued for the 



existence of a discrete hierarchy in fracture and rock properties |^6[, with a prefered 



scaling ratio around 3.5. Borodich ||6^ discusses the use of parametric-homogeneous 
functions for a parcimonious mathematical representation of structure presenting a 
hierarchy of log-periodicities. 



In their theoretical analysis. Ball and Blumenfeld predicted logarithmic os 



dilations in quasi-static crack growth, probably one of the very first example of such 
oscillations in non-tree structures. First, they coarse-grained a quasi-static growing 
crack as a wedge and found the behavior of the stress field around it. Then, they 
showed through a linear stability analysis that there is an instabitily to growth of 
branches where the stress is locally high. In this, it is conceptually similar to the DLA 
instability with respect to a locally high gradient of the field near the interface. Then 
they argued that because of this instability the dynamics couple to sub-dominant 
terms in the stress field to generate branches. So looking into what are the strongest 
subdominant terms, they could identify the behaviour of the branching. These terms 
happened to have the functional form 



On. — ^n, ' 



where an{< 0), and 6„ are numbers and Sn is the stress due to the contribution 
of the n-th term in the expansion of the stress field. Their argument was that, since 
these terms exhibit logarithmic oscillations, then so will the branching in the cracking 
pattern, resulting in a growing dendrite with logarithmically periodic branches. 

Let us finally mention that the distribution of avalanche strengths in the train 
model, a variant of the Block-Spring Burridge-Knopoff model of earthquakes pulled 
at one extremety, is a powerlaw times a log-periodic function |]70|. This reflects the 



spontaneous formation of a discrete hierarchy. In this model, the prefered scaling ratio 
in the length size of the avalanches is equal to 2. The moving blocks obey a discrete 
wave equation with attenuation at large velocities. A non-local cellular automata 
model retrieves the results and is defined so as to mimick the behavior of the train 
model. The main point is that the active sites of the cellular automata obey a discrete 
Laplace equation. In the overdamped approximation valid for large attenuation, the 



train model is also in this regime. The discrete hierarchy of the sliding events reflects 
the instabilities of the growth processes controlled by the Laplace equation 



5.3.2 Log-periodicity in time 



A growing body of evidences indicate that the log-periodic oscillations appear in 
the time dependence of the energy release on the approach of impending rupture in 
laboratory experiments ||22|, numerical simulations and earthquakes [|l^, |l^, 



It is thought that a similar type of cascade, from progressive damage 
at small scale to coalescence and unstable growth, is controlling the appearence of 
log-periodicity. The typical time-to-failure formula used in these works is 



Er^iu- ty 



1 + Ccos (^27r 



log(t^ - 1) 
log A 



(50) 



where E is the energy released or some other variable quantifying the on-going dam- 
age, tr is the time of rupture, m is a critical exponent, and \1/ is a phase in the 
cosine that can be get rid of by a change of time units. It has been found that 
the log-periodic oscillations enable a much better reliability of the prediction due to 
"lock-in" of the flt on the oscillating structure. Physically, the oscillations contain 
information on tr and thus help signiflcantly in its determination. A link between 
log-periodicity in space and in time is given in |7^. A typical fit by expression ( pD| ) 
to acoustic emission data is presented in figure 5. 



5.4 Cascade of structure in hydrodynamics 



Moffatt has studies similarity solutions for the fiow of a viscous fiuid near a sharp 
corner between two plates on which a variety of boundary conditions are imposed. 
For this, one has to solve the biharmonic Stokes equation for the stream function 
which admits separable solutions in plane polar coordinates {r,6): \|/ = fa{(^)- 
For angles between the two plates less than 146 degrees, the exponent a has been 
shown to be necessarily complex |j76[. Moffatt has shown that this result can be 



interpreted as implying the existence of an infinite sequence of eddies near the corner. 
It is interesting that viscosity, usually a damping mechanism, is here responsible for 
the generation of a geometrical progression of eddies. The damping has the effect, 
however, to give a large ratio (typically greater than 300) of the intensities between 
successive eddies. 

Recent experiments and analysis both numerical and analytical of droplet fission 
shows the existence of iterated instabilities that develop a discrete scale invariance. 
The reason for self-similarity is that, near breakoff, the droplet radius becomes much 
smaller than any other length scale, so that the shape of the interface becomes inde- 
pendent of these scales. Numerical simulations of the corresponding hydrodynamic 
equations with a weak noise source show that necks and blobs form repeatedly on 



smaller and smaller scales as the interface breaks |T3]. These instability cascades are 
only observed for ffuids with a viscosity greater than 1 P. The mechanism behind 
this cascade of instabilities is that, immediately before a neck forms, the thinnest sec- 
tion of the interface is well approximated by the similarity solution. The nonsteady 
singularity results from repeated instabilities of the similarity solution due to thermal 
capillary waves. 



5.5 Deterministic dynamical systems 



5.5.1 Cascades of sub-harmonic bifurcations in the transition to chaos 

An area where log-periodic structures should be expected is low-dimensional dynami- 
cal systems exhibiting the Feigenbaum sequence of subharmonic bifurcations to chaos 



|]78| . Indeed, this route to chaos can be understood from an asymptotically exact dis- 
crete renormalization group with a universal scaling factor. The existence of this 
prefered scaling ratio should thus lead to complex exponents and log-periodic oscilla- 
tions around the main scaling as the dynamics converges to the invariant Cantor set 
measure at criticality. It was noticed quite early that the length of the stable pe- 
riod diverges as a power law with log-periodic modulations as the control parameter 



approaches the transition to chaos. Argoul et al. |^ have studied the transitions to 
chaos in the presence of an external periodic field and show figures exhibiting very 
clearly that the Lyapunov exponent has a power dependence with log-periodic oscil- 
lations as a function of the amplitude of the external field. Similarly, the topological 
entropy at the onset of pruning in generalized Baker transformations is a power law 



function of the distance to the onset of pruning with log-periodic oscillations |8T 
The oscillations are due to the self-similar structure of the Cantor set forming the 
attractor |7^. 

5.5.2 Two-coupled anharmonic oscillators 

Consider the Hamiltonian 

H = ^ + ^{q'- br + ^ + ln^u^ + lgu {q^ - h^) , (51) 

where (P, u) and (p, q) are the two conjugate pairs of (momentum, displacement) for 
the two degrees of freedom u and q. The Hamilton- Jacobi equation of motion are 

'^ + aq{q^-h'')+guq = (52) 

+ = 0- (53) 

Suppose that q is observable while u is not. One can see the problem as the motion 
of a single particle in the time dependent potential 

KfT = ^{q' - br + u' + \9U {q' - b') , (54) 

where u{t) is a time-dependent parameter responding to the particle position q{t). 
The quartic potention Ves exhibits two wells and the particle oscillates between in 
a chaotic fashion for sufficiently large coupling g. West and Fan have shown 



that this system is chaotic with positive Lyapunov exponent for sufficiently large 
coupling g. More interestingly, the spectrum of q{t), equal to the Fourier transform 
of the correlation function {q(t)q{t + r)), is found to be a power law ~ with log- 
periodique modulations. These modulations are best seen in the so-called "rescaled 



range analysis" introduced by Hurst |^ , which emphasizes the extremal ffuctuations 
of the signal and thus the observation of the modulations. The log-periodic oscilla- 
tions disappear progressively at large time scales, due probably to the fact that they 



constitute a correction to scaling. Intuitively, one can propose the following qualita- 
tive understanding. The time-dependence of u makes the dynamics chaotic and also 
modulates the height of the barrier separating the two wells close to ±6. For a given 
height and in the presence of chaotic behavior, one can expect a kind of Arrhenius 
activated behavior controlling the time of passage from one well to the other. The 
log-periodicity would then result from a discrete possible sets of values of u at which 
this cross-over occurs. Further analysis is needed to ascertain this mechanism. 



5.5.3 Near-separatrix Hamiltonian chaotic dynamics 



In a series of papers, Zaslavsky et al. (see for a review) have introduced a renor- 
malization group transform for Hamiltonian dynamical systems with a saddle-point 
(separatrix) presenting a logarithm dependence of the period of periodic orbits as a 
function of the distance in energy scale to the saddle-point. The scale transforma- 
tion performed in the renormalization group analysis corresponds to change both the 
strength e of the perturbation applied to the system and the distance h to the saddle- 
point. In analogy with usual critical points, e plays the role of the scale L at which 
the systems is coarse-grained and h is the distance to the critical point. For the har- 
monically perturbed pendulum and for the motion of two particles in a fourth-order 
polynomial potential, the renormalization group involves a prefered scaling ratio A 



for the scale transformation (eq.3.17 and eq.4.14 in |p4[)- the language of this re- 
view, this corresponds to having a discrete scale invariance. The consequence is that 
suitable observables should have complex critical exponents and the signals should 
be decorated by log-periodic signatures. 



5.6 Animals 

We have notice 0] that, in contrast to common lore, complex critical exponents 



should generally be expected in the field theories that describe geometrical systems, 
because the latter are non unitary. In particular, evidence of complex exponents in 
lattice animals, a simple geometrical generalization of percolation has been presented 



T0| . The model of lattice animals is the most natural generalization of the percola- 
tion model, which itself is the prototype of disordered systems. The animal problem 
is the statistics of connected clusters on a lattice |^ ^ and thus also describes un- 



rooted branched polymers. Using transfer matrix techniques, the number of unrooted 
branched polymers of size N is found to exhibit a correction to the main scaling with 
a complex exponent. 

Recall that in the percolation problem, bonds are occupied with a probability p 
and unoccupied with probability 1—p. For a given configuration, connected parts are 
called clusters. To study the statistics of one percolation cluster, one can sum over all 
configurations for the bonds that do not belong to this cluster nor to its perimeter. 
Since they are either occupied or inoccupied, the sum over all configurations just gives 
a unit weight. Hence in percolation, clusters are simply weighed with — p)^^ 

where A^f, is the number of bonds in the cluster, Np the number of bonds of the 
perimeter. Now the animal problem is a more general model where a cluster is 
weighed by p^^qNp -yyii;}^ general values of p, q. By varying p,q a critical point is 
met which is always in the same universality class of so-called animals. Only when 
p + q = 1 is this critical point in a different universality class, percolation, which 



therefore can be considered as a tricritical point in the animals parameter space. 

The result of the analysis of |0 is that the number T/v of unrooted branched 
polymers of size N in the plane is given by 



N 



^K2-2Xi)-3 ^ ^^^(2-2X2,H-3) (2X2,!^ \og N + 



(55) 



where u is the radius of gyration exponent, u ^ 0.64. Recall that the leading term 
in (|^) is actually known exactly z/(2 — 2Xi) — 3 = — 1 [^. Hence we see log- 
periodic terms to appear in the next to leading behaviour of T/y. Unfortunately, 
since conformal invariance is broken, this argument does not allow us to make any 
predictions on the amplitude of these terms, which might well be very small. The 
DLA problem is much more favorable in that respect probably due to enhancement 
effects stemming from the long-range interactions of Laplacian fields. 



5.7 Quenched disordered systems 

Renormalization group analysis of a variety of spin problems with long-ranged quenched 
interactions have found complex critical exponents ISr 



89, 90, 91, 92 



However, these 

authors have in general remained shy as to the reality of their results. Indeed, one 
could argue that uncontrolled approximations (present in all these works) rather than 
physics could be the cause of the complex exponents. Derrida and Hilhorst have also 
found log-periodic corrections to the critical behavior of ID random field Ising model 



at low temperature by analyzing products of random non-commutative matrices ||93|| . 

With the qualitative understanding of the ultrametric structure of the energy 
landscape of spin glasses [0, ^ in the mean field approximation, one could conjecture 
that these above results could be the observable signature of the hierarchical structure 
of energy states in frozen random systems. The problem is that, even if hierarchical, 
the ultrametric structure is believed to be continuous and it is not clear what could 
produce the discrete scale symmetry. It is generally believed that such topology occurs 
more generally in other complex systems with highly degenerate, locally stable states 
96|. However, much works remain to be done to clarify this problem. The additional 



presence of long-range interactions complicate the matter further. 

A dynamical model describing transitions between states in a hierarchical system 
of barriers modelling the energy landscape in the phase space of meanfield spinglasses 
leads again to log-periodic corrections to the main logt behavior [1^. 

Let us also mention that mth critical Ising models (m = 1 for the Gaussian model, 
m = 2 for Ising, m = 3 for tricritical Ising, ...) have a free energy exhibiting log- 
periodic oscillations as a function of the control parameter for large m, a signature 
of the geometrical cascade of multicritical points P3, 0^, [10 . 



6 OTHER SYSTEMS 

6.1 The bronchial tree 

It has been pointed out that the morphology of the bronchial airway of the mam- 
malian lung is roughly hierarchical leading to a log-log plot of the average diameter 
of a branch of a mammalian lung (for human, dog, rat and hamster) as a function 



of the branch order which exhibits a full S-oscillation (log-periodic) decorating an 
average linear (power law) dependence. This fractal structure has been argued to 
allow the organ to be more stable with respect to disturbance ||100| , |101| , |102| , |10^ 



but the physical mechanism underlying its appearence is not understood. 
6.2 Turbulence 

Probably the first theoretical suggestion of the relevance of log-periodic oscillations to 
physics has been put forward by Novikov to describe the influence of intermittency in 
turbulent flows ||^. The idea is that the DSI could stem from the existence of a pref- 
ered ratio in the cascade from large eddies to small ones. The existence of log-periodic 
oscillations has not been convincingly demonstrated as they seem quite elusive and 
sensitive to the global geometry of the flow and recirculation |p.04] , |105|| . Shell models 



of turbulence, which have attracted recently a lot of interest pO(j[j construct explicitely 



a discrete scale invariant set of equations. In these models, self-similar solutions of 
the cascade of the velocity fleld and energy in the discrete log-space scale have been 



unravelled ||107|| , whose scaling can be related to the intermittent corrections to Kol- 
mogorov scaling. We note that some of these solutions rely on the discrete scaling 
shell structure and would disappear in the continuous limit. However, the relevance 
of these discrete hierarchical models and more generally of log-periodic oscillations 
have not been explored systematically and their conflrmation in turbulence remains 
open. 

Based on theoretical argument and experimental evidence, the author has con- 
jectured that structure functions of turbulent times series exhibit log-periodic mod- 



ulations decorating their power law dependence [108|. In order to provide ironclad 
experimental evidence, we stress the need for novel methods of averaging and propose 
a a "canonical" averaging scheme, amounting to condition the average of structure 
factors of turbulent flows on local maxima of the energy dissipation rate or the energy 
gradient. The strategy is to determine the scale at which the dissipation rate is 
the largest in a given turn-over time series. This speciflc scale Vc translates into a 
speciflc "phase" in the logarithm of the scale which, when used as the origin, allows 
one to phase up the different measurements of a structure factor Sp{r) = Ap{ery^^ in 
different turn-over time realizations. We expect, as in Laplacian growth and in rup- 
ture, that the log-periodic oscillations will be reinforced by this canonical averaging. 
Demonstrating unambiguously the presence of log-periodicity and thus of discrete 
scale invariance (DSI) in turbulent time-series would provide an important step to- 
wards a direct demonstration of the Kolmogorov cascade or at least of its hierarchical 
imprint. 

6.3 Titius-Bode law 



Dubrulle and Graner ||109|| have noticed that the Titius-Bode law of planets distance 



to the sun r„ = roK^ with K ^ 1.7 can be seen a discrete scale invariant law {K 
then plays the role of A in our notation). They show that all models that have 
been proposed to explain the Titius- Bode law share the common ingredient of scale 
symmetry. Assuming a discrete symmetry breaking in the rotation invariance, they 
thus show that any such mechanism is compatible with the Titius-Bode law. As a 
consequence, this law cannot a priori be used to constraint the mechanism of planet 



formation and their organization around the sun. What is however not understood 
is the physical mechanism, if any, at the basis of the breakdown of continuous to 
discrete scale invariance embodied in the Titius-Bode law. 



6.4 Gravitational collapse and black hole formation 



Choptuik [|110|| has recently shown that, in contrast to the general view, black holes of 



mass smaller than the Chandrasekkar limit could be formed and that, in the process of 
formation, the solutions would oscillate periodically in the logarithm of the difference 
between time and time of the formation of the singularity. This gravitational collapse 
is an example of critical behavior, describing how the mass M of the black hole 
depends on the strength p of the initial conditions: M {p — p*)^ for p > p* and 
otherwise, where p* is the threshold value. It has been shown that classical close- 
to-critical black holes (obeying Einstein's equations) coupled to a massless complex 
scalar field have a leading real exponent and a sub leading complex exponent [ 11CI|| , 



which would correspond to a log-periodic spectrum of masses. Alternatively, the real 
and complex exponents control the time development of the black hole instability 
which is also log-periodic in time, corresponding to continuous phase oscillations of 
the field. 



6.5 Spinodal decomposition of binary mixtures in uniform 
shear flow 



Corberi et al. have studied the phase separation kinetics of a binary mixture 



subjected to an uniform shear flow quenched from a disordered to a homogeneous 
ordered phase. They used the time-dependent Ginzburg-Landau equation with an 
external velocity term. They flnd numerically that the typical sizes of domains in 
the directions parallel and orthogonal to the direction of the shear as well as the 
excess viscosity exhibit log-periodic oscillations as a function of time that decorate 
the scaling laws. They have discovered this behavior from the numerical integration 
of the Langevin equation for the evolution of the system. From a theoretical point 
of view, the authors do not have an explanation for the log-periodic behavior of 
the scaling functions (private communication). They interpret these log-periodic 
oscillations as due to a growth mechanism where stretching and break-up of domains 
occurs cyclically in In t. They have performed a one-loop perturbation approximation, 
which together with the fact that (1) the natural time variable is Int and (2) the 
scaling function is assumed periodic allows them to retrieve the log-periodicity of their 
numerical simulations. This is a flrst step, but it does not provide a mechanism for 
log-periodicity since it has been "put by hand" in the scaling ansatz. It is particularly 
interesting that the alternative stretching and break-up lead to this cascade of scales. 
It would be of interest to pinpoint the underlying physical mechanism by while a 
prefered scale ratio is selected. 

6.6 Cosmic lacunarity 

The present observable distribution of galaxies in space provides a signature of their 
formation and interaction. There is an undergoing controversy on the nature of 
the spatial distribution of galaxies and galaxy clusters in the universe with recent 



suggestions that scale invariance could apply ||112|| . Recently, evidence of "cosmic 



lacunarity" has been presented ||113|| : the prefactor of the power law dependence 
of the gth moment of the spatial galaxy distribution oscillates periodically with the 
logarithm of the distance r with the same periodicity for g = 2, 3 and 4. The corre- 
sponding scaling factor is around 2^'^ ^ 3. A similar log-periodic oscillation has been 
suggested much earlier ||114|| with approximately the same scaling ratio. F. Combes 



(private communication) cautions that all previously announced periodicities (not 
log-periodicity) have finally been proven wrong. There are many bias in the data and 
we have access only to the luminosity, not directly to the masses. Furthermore, F. 
Combes argues that, since the existence of a prefered scaling ratio needs to "nucle- 
ate" from either a small scale (like in DLA where the small scale is the size of the 
particles) or from a large scale (maybe like turbulence), it is not obvious what should 
be these cut-off scales in the present case. 

6.7 Rate of escape from stable attractors 

Let us also mention the recently discovered log-periodic behavior of the rate of escape 
from a stable attractor surrounded by an instable limit cycle as a function of the 



strength of the white noise |115|. This is an example where the rate of escape, as 



calculated from a Fokker-Planck equation, is non-Arhenius. 
6.8 Interface crack tip stress singularity 

Complex singularities are also found in the divergence of the stress as a function of 
the distance to the tip of a crack at the interface between two different elastic media 



116|| . The standard a ~ r 2 singularity, where r is the distance to the crack tip and 
a is the stress, is replaced by 

a ~ r-^+^'^ , (56) 



where 

1 



Subscript 1 and 2 refer to the material in y > and y < 0, respectively; k = 3 — Au 
for plane strain and k = |^ for plane stress, u is the Poisson ratio and /i is the 
shear modulus. Interface cracks have important practical applications since inter- 
faces between composite media are often the locii of damage nucleation leading to 
the incipient rupture. The existence of this complex singularity suggests that the 
mechanism of damage and rupture at interfaces could be quite different from that 



in the bulk ||117|| . Following this work, a wealth of studies have followed (see [|117 
and | 118 | and references therein), but the physical understanding of the appearance 
of a complex critical singularity has remained elusive. Since the solution shows that 
the two modes of deformations in tension and shear (modes I and II) are intrinsi- 
cally coupled for an interface crack in constrast to what happens for a crack in an 
homogeneous medium, one could hope to identify the physical origin of the com- 
plex exponent in this coupling. Let us also mention that the pressure distribution 
as a function of distance to the corner in the Hertz problem of two different elastic 
spheres compressed against each other is also described by a power law with complex 
exponent ||119|| . 



See also section 5.3.1 on Log-periodicity in the geometrical properties, which refers 
to the work of Ball and Blumenfeld . 



6.9 The Altes'wavelet and optimal time-scale product 

The Altes wavelet has a log-periodic structure, with a local frequency varying hyper- 
bolically. It has the remarkable property of minimizing the time-scale uncertainty, in 
the same sense as the Gaussian law minimizes the time-frequency uncertainty. As will 
become clear below, it also has the nice property that differentiating it corresponds 
to dilate it by a fixed factor. 

Consider a signal U{t) with Fourier transform U{uj). Let us define 

1 f+oo 

En = — du liju'r U{uj)\^ (58) 



27r J- 



oo 



as the energy of the nth signal component. Then 



Vn{uj) = 7= (59) 



is an energy-normalized signal component. The rational to study Vn{uj) is that it 
corresponds to the normalized spectrum of the nth derivative of the signal U{t). The 
analysis of the signal U{t) through the different n components Vn{uj) corresponds to 
analyse the different derivatives of the signal. 

Let the mean-square bandwidth of the function (ju)"- U{uj) be 

" !l^du,\{juj,YU{uj,)\^ ■ ^ ' 

It is easy to verify from the Schwarz inequality that [|120| 



> Wl, . (61) 

This result means that the band-width increases upon differentiating the function 
U{t). Altes proposed to measure signal and filter complexity in terms of time- 
bandwidth product T W. Equation ( |6TD says that increasing n will result in a larger 
W, and in many cases T W will also increase. The components Ki(a;) will thus 
generally become more complicated and more difficult to analyse as n increases. 
The idea of Altes ||1 20|| is to postulate a signal with constant T W : 

T„ PV„ = T„_i Wn-i . (62) 

This signal must have a root-mean-square duration T that decreases when the time 
function is differentiated. The solution is that Vn{uj) satisfies 

Vn{uj) oc Vn-i{uj/k) , wherc k > 1 . (63) 



Indeed, from (|63|) , Wn = k Wn-i and T„ = T„_i//c. The complexity of the compo- 
nents is the same : the inevitable bandwidth increase is compensated by a propor- 
tionate decrease in duration. 



It follows from (^) that 



(64) 



leading to 

cj" Uiu) OC f/(cj/r) . 
Specifying the constant which can depend on n leads to 

cj" Uiu) = C(n) [/(o^/r) . 



(65) 



(66) 



This functional equation is very similar to the renormalization group equation previ- 
ously encountered. The difference stems essentially from the fact that the constant 
C{n) can depend on n. In the present approach, instead of changing the scale, we 
change the order n of derivation of the signal, which corresponds to changing the 
weight of the components. This corresponds to putting more and more weight to 
the high frequencies, and thus, decreasing n is like a coarse-graining. Then, the con- 
dition (|6^) of equal complexity for all components provides the self-similar function 
equation (p6|). The solution of this functional equation is ||120 



U (uj) = A uj" exp 



(lncj)2 
2 In A; 



/ 27rc \nui' 



where 



A is a constant and 



u = Re{^-^}, 271 c = Im{C'{0)} 
k 



C{n) = . 



(67) 

(68) 
(69) 



The origin of the log-periodicity of the waveform given by (|67|) is clear : it is the con- 
junction of the fact that derivative orders n are integers and that each additional order 
of the differentiation acts exponentially on the frequency spectrum as proportional to 
uj"'. This example is thus again a conversion of discrete translation to discrete scaling 
by an exponentiation mechanism. 

U (t) has a waveform which is a function of the parameters i>,k,c. The parameter 
k determines the spectral width, a larger k implies a greater bandwidth, k also 
determines how much stretched is caused when the Fourier transform is multiplied 
by [juj)^, i.e. under n derivatives of U{t). The parameter c determines the slope of 
the spectral phase function and is related to the "chirp" that appears in the times 
waveform. A larger value of c leads to more zero crossings. It can be shown that 
the waveform (^) is optimum as having a lack of Doppler coupling and a maximum 
Doppler tolerance [|120|| . Note that, by construction, all the components Vn{uj) have 
the same ratio of center frequency to bandwidth. Remarkably, these waveforms ([67| ) 
are very similar to those employed by animal sonars, such as bats and dolphins. 



6.10 



Orthonormal Compactly Supported Wavelets with Op- 
timal Sobolev Regularity 



Ojanen ||121|| has recently presented a new construction of an orthonormal compactly 
supported wavelets with Sobolev regularity exponent as high as possible among those 



mother wavelets with a fixed support length and a fixed number of vanishing mo- 
ments. The increased regularity is obtained by optimizing the locations of the roots 
the scaling filter has on the interval [7r/2,7r]. The obtained wavelets have a discrete 
scale invariant structure, as is apparent from their examination. This is also seen 
directly by their construction. 

To construct orthonormal compactly supported wavelets, Ojanen uses as the start- 
ing point the discrete dilation equation 

= V2J2 Ck(t>{'^x -k), X eR, CkE R, (70) 

kez 

for the scaling function or father function 0. A solution exists in the sense of distri- 
butions when I^fcCfc = and it has compact support if only finitely many of the 
filter coefficients are non-zero (in fact, if Cfc 7^ only for k G {0, . . . ,n} then the 
support of (f) is contained in [0, n]). The solution of the dilation equation is given by 

00 

0(0 = n ^0(2-^0, (71) 

where the trigonometric polynomial mo(C) is the scaling filter of (p, 

^oiO = ^ E c.e-^'^^. (72) 



The condition when the solution to ^ yields wavelets can be expressed in terms 
of the scaling filter mo(0- If ^0 satisfies 

l^o(Or+l^o(vr + OI' = l, eei?, (73) 

mo(0) = 1, and a technical condition called the Cohen criterion holds, then there 
exists a scaling function and wavelet pair. 

If is such that it gives rise to a wavelet, the wavelet ip is obtained by 

V^(a^) = E(-l)'ci-fc<^(2a;-^)- (74) 

kez 

Note that then and ip have the same regularity properties. 

The original construction of orthonormal compactly supported wavelets by Daubechies 
||122|| considers a trigonometric polynomial rriQ with at most 2N non-zero coefficients 
Co, . . . , C27V-1- The polynomial mo is required to satisfy |7^ and to have a zero at 
vr of maximal order (the order of the zero is also the number of vanishing mo- 
ments the corresponding wavelet has). The resulting equations can be solved and the 
constructed mo has no roots on (— 7r,7r), hence it satisfies the Cohen criterion. 

To construct smoother wavelets, Ojanen's approach is to reduce the order of the 
zero mo has at vr and instead to introduce zeros on the interval [it/ 2, it]. Keeping 

and the number of zeros on [7r/2,7r] fixed, Ojanen uses numerical optimization 
to choose the locations of the roots so that the resulting wavelets are as smooth as 
possible. He is thus able to construct wavelets that are more regular than those 



introduced in 122, 123, 1241. 



6.11 Eigenfunctions of the Laplace transform 



Log-periodicity and complex exponents play a very important role in integral equa- 
tions of the type (?(r) = /q°° K{vT)p{v)dv, with < r < oo where the kernel K has 
the property \K(x)\x~^^'^dx < oo. This class of equation includes the Laplace 
transform, the Fourier sine and cosine transforms and many other integral equa- 
tions of importance in physics. It is notorious that the inversion problem of getting 
p{v) from the measurement of g{T) is ill-conditioned. This can be seen to result 
from the form of the eigenfunctions and eigenvalues of the Laplace transform and 
similar dilationally invariant Fredholm integral equation [|126|| . For instance, the 
eigenfunctions of the Laplace transform, which form a complete orthogonal basis, are 
(p^iv) = v~^^'^ cos{ujlnv — 9aj) and = — f ^^^^ sin(c<jlnt> — 9^), where 6*^ is a 

function of u. The eigenvalues are exponentially decreasing with u and this controls 
the ill-conditioned nature of the Laplace inversion. The log-periodicity of the eigen- 
functions lead to an optimal sampling determined by a generalized Shannon theorem 
which obeys a geometrical series ||127| . 



6.12 Life evolution 



A very intringuing study has been performed by Chaline et al ||128|| , who have found 
that the fossil equine of North America, the primates and the rodents have followed 
an evolution path puntuated by major events that seem to follow a geometrical time 
series of the type (|75D, with a prefered time scale factor A ~ 1.7. Using t-student 
statistics and Monte-Carlo tests, Chaline et al. find a very high statistical significance. 
The critical time tr is roughly the present for the equines, in agreement with the 
extinction of this species in North America ten thousand years ago (which is very 
short compared to geological time scales), about two millions years in the future for 
the primates and between 10 to 60 millions years in the future for the rodents. The 
main problem is to ascertain that the determination of the times of the major events 
has been done independently from the mathematical analysis. Otherwise, this leads 
to a well-known statistical bias and the statistical significance of these results looses 
its meaning. 



7 APPLICATIONS 

7.1 Identifying characteristic scales 

In our opinion, the main interest in identifying log-periodicity in data is the charac- 
terization of the characteristic scales associated to it. Indeed, it must be clear that 
the log-periodic corrections to scaling imply the existence of a hierarchy of character- 
istic scales (in space or time). For instance, in the time-to-failure analysis given by 



(|50|), the hierarchy of time scales is determined by the local positive maxima of the 



function E. They are given by 

tr - t„ = rAt , (75) 

where r oc exp (-^tan-i^;^). The spacing between successive values of t„ ap- 
proaches zero as n becomes large and t„ converges to tr- This hierarchy of scales 
tr — tn are not universal but depend upon the specific geometry and structure of the 



system. What is expected to be universal are the ratios ^^^^ = A 2 . From three 
successive observed values of t„, say tn, tn+i and t„+2, we have 

tr = ^'+1 ~ ^"+2^" . (76) 
2tn+l ~ tn — tn+2 

This relation applies the Shanks transformation to acceleration of convergence of 
series. In the case of an exact geometrical series, three terms are enough to converge 
exactly to the asymptotic value t^. Notice that this relation is invariant with respect 
to an arbitrary translation in time. In addition, the next time t^+s is predicted from 
the first three ones by 



These relations have been used in |35, 21, UW. Physically, time or space scales give 



us access to additional information and clues about the underlying processes and the 
existence of a hierarchy of prefered scales, as in DSI, will tell us something about the 
underlying processes. This is lost in usual critical behavior in which all scales are 
treated as playing the same role. 

7.2 Time-to-failure analysis 

Another important application of log-periodicity is its use in making more robust 
and precise time-to-failure analysis. We have already mentionned the importance of 



log-periodicity for predictions |18|, |19|, |21|, The derived time-to- failure 



analysis is now being implemented for routine industrial testing in the space industry 
in Europe. As already mentionned, t he reason for this improvement is that a fit can 
"lock-in" on the oscillations which contain the information on the time of failure and 
thus lead to a better prediction. 

7.3 Log-periodic antennas 

Let us mention the engineering application of antennas using log-periodic electro- 
magnetic antennas ||129| , |130| , |131|, |132| , [133| , |134| , |135|| . The DSI structure provides an 



optimal compromise between maximizing bandwidth and radiation efficiency. 
7.4 Optical waveguides 

Graded-index optical waveguides with optimized index profiles can support a family 
of weakly localized modes with algebraic tails with log-periodic modulations in the 



evanescent field ||125|| . The log-periodic oscillations result from an interplay between 
the critical nature of the modes and absorption (complex index of refraction). This 
could have applications in techniques using evanescent waves. 



8 OPEN PROBLEMS 



8.1 Nonlinear map and multicriticality 



In this brief review, we have kept the analysis of the renormahzation group at the level 
of a linear expansion of the flow map. Taking into account the nonlinear structure 
of the flow map, as for instance in (^), may lead to an infinite set of singularities 
accumulating at the main critical point or even to the whole axis of the control 



parameter being critical (in the chaos regime of the flow map [p,36|| ). The question of 
the relevance of these regimes to nature is still open (see |]19[ for a proposed application 
to earthquakes). Generalization to several control parameters and multicritical points 
would be useful. 



8.2 Multilacunarity and quasi-log-periodicity 



We have seen that DSI embodies the concept of lacunarity. The set of complex expo- 



nents or singularities has been called multilacunarity spectrum 1^^. Generalizations 
with several different incommensurate log-frequencies would be of great interest and 
seem to appear for instance in the DLA problem |T^. Complex multifractal dimen- 
sion spectrum in the presence of disorder can be handled using probabilistic versions 
of the renormahzation group [ |137[ ] and their development and impact are just emerg- 
ing. Note also that the set of complex exponents provides a better characterization 
of the underlying multiplicative process and could improve the conditionning of the 
inverse fractal problem ||138|| . The g-derivative is a natural tool to discuss homoge- 
neous functions with oscillatory amplitudes. It has recently been used to describe 
cascade and multifractal models with continuous scale changes [p.39|| . In Ref. |p.4q] , a 
kinetics is built upon the g-calculus of discrete dilatations and is shown to describe 
diffusion on a hierarchical lattice. 



8.3 Effect of disorder 

A very important practical question is the effect of disorder and the process of aver- 
aging. Disorder is expected to scramble the phases of the log-periodic oscillations and 
it is a priori not clear whether the log-periodic oscillations are robust. It turns out 
that small fluctuations around the log-periodic structure do not seem to spoil DSI, 
as found in many examples quoted above. For instance, in the needle DLA problem, 
intervals between needles were taken to fluctuate by a few percent without altering 
signiflcantly the log-periodic structure of the growth process |^0[ . DLA clusters them- 



selves are formed under a very strong annealed noise, corresponding to the random 
walk motion of the sticking particles. Nevertheless, clear evidence of log-periodicity 



in the mass as a function of radius has been found |T5[ giving confldence in their 
robustness with respect to disorder. This has been further substantiated by explicit 
calculations |]10[ showing generally that the complex exponents are robust. 



8.4 Averaging : grand canonical versus canonical 

However, we must stress that disorder introduces a sensitive dependence of the phase 
in the cos log formula : different realizations have a different phase and averaging will 
produce a "destructive interference" that makes vanish the log-periodic oscillations. 
It is thus important to carry out analysis on each sample realization separately, 
without averaging. For instance, in the DLA case, 350 clusters of 10^ particles have 



each been analyzed one by one and an histogram of the main log-frequencies has been 
constructed. Theoretically, preventing averaging is a problem as one is usually able 
only to calculate quantities averaged over the different realizations of the disorder. 
However, it must be stressed that the fact that log-periodic oscillations are mainly 
present before averaging tells us that they are specific fingerprints of the specific 
system one is looking at. This is obviously a desirable property for prediction purposes 
in engineering and other practical applications. An open problem however is to devise 
optimal tools to decipher the log-periodic structures in highly noisy data, as is usually 
the case due to limitation of sizes for instance. We note also that going to very large 
systems will in general progressively destroy the log-periodic structures as they are 
often correction to scaling |. Random versions of Cantor fractal sets have recently 
been shown to exhibit robust log-periodic structures even when averaging ||141| . 



Pazmandi et al. [|144|| have recently argued that the standard method of averaging 
carried out in disordered systems introduces a spurious noise of relative amplitude 
proportional to the inverse square root of the system size. This so-called "grand 
canonical" averaging can thus destroy more subtle fluctuations in finite systems, 
controlled by a correlation length exponent less than | (| is the minimum value of 
the correlation length exponent that would not be hidden by the usual grand canonical 
averaging). Pazmandi et al. thus propose an alternative averaging procedure, the 
so-called "canonical" averaging, which consists in identifying, for each realization, 
the corresponding specific value of the critical control parameter Kj^. The natural 
control parameter then becomes A = {K — K^)/K^ and the act of averaging can 
then be performed for the samples with the same A. We have used this procedure to 
identify log-periodicity in the elastic energy E prior to rupture in a dynamical model 
of rupture in heterogeneous media. Previous works have shown that E follows a power 
law E {tr — as a function of the time to failure ||145|| . Performing the usual 



(grand canonical) averaging over twenty different realizations of the disorder provides 
very good evidence of the power law but no evidence of log-periodicity. We have thus 
developed the following alternative averaging procedure ||146|| . We constructed the 
second derivative of E with respect to time for each realization, thinking of it as 
a kind of susceptibility. The time at which this second derivative is maximum 
has been identified and this point has been used as the effective value of the time tr 
of rupture for each realization. Then, the first derivative, giving the rate of energy 
released, is averaged over all samples with the same (t^—t)/t^. The result is presented 
in figure 6 in a log-log scale. Four to five approximately equidistant spikes (in log 
scale) are clearly visible. The log-period allows us to identify a prefered scaling ratio 
A = 2.5 ± 0.3. It is probable that similar averaging procedures better tailored to get 
rid of spurious fluctuations from realization to realization will play an increasing role 
in the physics of disordered and turbulent media. 



8.5 Amplitude of log-periodicity 

Log-periodicity is found in spin systems in hierarchical lattices. However, the effect is 
usually very small, typically 10~^ or less in relative amplitude. In contrast, we have 
found it much stronger in rupture and growth process, typically 10~^ or so in relative 



9 



This is often due to the disorder which can be shown to rcnormahzc the real part of the complex 



exponents so that they correspond to sub-leading correction to the main scaling behavior 10, 



amplitude. The reason is not very well understood but might stem from the strong 
amplification effects occurring in such Laplacian fields. 

There is a situation where the factors controlling the log-periodic amplitudes are 
understood. This is the case for spin models in hierarchical networks. In this case, 
the free energy of the high temperature phase is in general analytic close to in a 
non-vanishing angular sector in the complex plane. As a consequence, the Fourier 
coefficients Cn of the log-periodic amplitudes decay exponentially with their order 
n. In practice, the first coefficient Ci is already given by the exponential law and 
therefore very small (see eq.(28) in In many other cases, the coefficients c„ 

decay more slowly, typically as powers of n, such that the Fourier sum which controls 
the amphtude of the log-periodicity does not possess a domain of analyticity. This 
occurs for instance for the density of states of harmonic chains with random masses, 
both for the Lifshitz singularities at the band edges or at the Halperin singularities 
(cluster resonances) in the case of binary mixtures ||142|| . There are other cases 



where the log-periodic amplitudes may be of both types depending on what we look 
at : in a chain of quantum XY spins with quasiperiodic modulated couplings, the 
singularities of the density of states correspond to the first type (exponential with 
domain of analyticity) while the singularities corresponding to the magnetization and 
the susceptibility are of the second type (power decay and no domain of analyticity). 
See [P3 compare figures 3 and 4 to figures 5 and 6. 



8.6 Where to look for log-periodicity 

From the point of view of non-unitary field theory, we should expect generically the 
existence of complex exponents. However, there is no known recipe to tell us what are 
the relevant observables that will have complex dimensions. Practically, this means 
that it is not a priori obvious what measure must be made to identify log-periodicity. 
In other words, it is important to look carefully at available data in all imaginable 
angles to extract the useful information. Of course, one must always been aware of 
statistical traps that noise can be taken for log-periodicity. Analysis must thus be 
carried out with synthetic tests for the null hypothesis, bootstrap approach, etc (see 
for instance |2T|, ^ for the application of statistical tests and the bootstrap method 
in this context). 



8.7 Prefered scaling ratio around 2? 

Another puzzling observation is the value of the prefered ratio A ~ 2, found for a 
wide variety of systems, such as in growth processes, rupture, earthquakes, financial 
crashes. H. Saleur (private communication) has noticed that 2 is in fact the mean field 
value of A obtained by taking an Ising or Potts model (with Q states) on a hierarchical 
lattice in the limit of an infinite number of neighbors. Consider a diamond lattice 
with n bonds connected to the upper and lower nodes (the usual diamond lattice 
discussed above has n = 2). The discrete renormahzation group equation connecting 
K = e^"^, where /? is the inverse temperature and J the coupling coefficient, from one 
generation to the next is the generalized version of (35) to n-bonds connectivity: 



K'(K) = lf(K)r = (^if^)" . 



In the limit n — cx) where the number of coupled nodes increases without bounds, 
we expect physically the ordered-disordered transition to occur at larger and larger 
temperature, corresponding to a fixed point of (^) K'{K*) = K* — >■ 1. Expanding 
around 1, we indeed find K* = 1 + — asymptotically. The linearization of the 
renormalization group map (^) gives K' — K* = X{K — K*) with A = nK* '^^^^^ \k* 
2 in the limit n ^ oo. Can this argument be extrapolated to out-of-equilibrium 
systems? 

8.8 Critical behavior and self-organized criticality 

Time-to-failure analysis of earthquakes seem at variance with the globally stationary 
view point, for instance captured by the concept of self-organized criticality as applied 
to plate tectonics ||143|| . Recently, it has been shown ||74| that a simple model of earth- 
quakes on a pre-existing hierarchical fault structure exhibits both self-organization 
at large times in a stationary state with a power law Gutenberg-Richter distribu- 
tion of earthquake sizes. In the same token, the largest fault carries irregular great 
earthquakes preceded by precursors developing over long time scales and followed by 
aftershocks obeying the j Omori's law of the rate of seismicity after a large earth- 
quake. The cumulative energy released by precursors follows a time-to-failure power 
law with log-periodic structures, qualifying a large event as an effective dynamical 
(depinning) critical point. Down the hierarchy, smaller earthquakes exhibit the same 
phenomenology, albeit with increasing irregularities. The study of the robustness of 
this scenario for other models and situations is an open question. 
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FIGURE CAPTIONS: 



Figure 1 : Construction of the triadic Cantor set : the discrete scale invariant 

geometrical structure is built by a recursive process in which the first step consists 

in dividing the unit interval into three equal intervals of length \ and in deleting the 

•J , 

central one. In the second step, the two remaining intervals of length ^ are themselves 
divided into three equal intervals of length g and their central intervals are deleted, 
thus keeping 4 intervals of length i, and so on. 



Figure 2: Construction of the hierarchical diamond lattice used in the Potts model. 
This lattice is obtained by starting with a bond at magnification 1, replacing this bond 
by four bonds arranged in the shape of a diamond at magnification 2, and so on. The 
spins are placed at the sites. At a given magnification 2^, one sees AP bonds, and thus 
|(2 + 4P) sites. 



Figure 3: 



where z/^j 



log(l-p) 



is the real part of (^) as a function of Int. 



t^^R ' - --ft log{u/v) 

The averaging has been performed over different realizations of the random walk 



(taken from [jSO 



Figure 4: a) Map of 5000 needles which have grown according to the DLA rules 
from an initial configuration where all the needles were approximately of the same 
length equal to their average separation. We have used a periodic lattice and added 
a small random value (a few percent of the period) for their lateral position. This 
configuration corresponds to the time when the largest needle has a length equal to 
one-third of the size of the system, b) The probability density function of the needle 
lengths shown in a) in where very clear log-periodic oscillations decorate the power 
law behavior, (see pO| ). 



Figure 5: Logarithm of the acoustic emission energy released as a function of the 
pressure (in bars) applied within a pressure tank made of matrix-fiber composite at 
the approach of rupture. The continuous lines correspond to the best fit by expression 
dg) (see |22|). 



Figure 6: Rate of elastic energy released as a function of the logarithm of the time 
to failure in the dynamical model of rupture with damage introduced in ||145|| . The 
dots have been obtained using the "canonical" averaging procedure discussed in the 
text. 



